Circulación de fluidos a través Q 
de lechos porosos. Fluidización M 


9.1. Introducción 

El flujo de fluidos a través de lechos de partículas sólidas aparece en numero- 
sos procesos técnicos. Así, en operaciones de filtración y flujo a través de colum- 
nas de relleno, utilizadas en operaciones de destilación, absorción, adsorción e in- 
tercambio iónico, se presentan ejemplos de este tipo de circulación. En el caso de 
la filtración, las partículas sólidas se depositan en un medio filtrante, y el fluido 
pasa a través de la masa filtrante depositada. En otros procesos, como la adsorción, 
intercambio iónico y otros, el fluido se mueve a través de un lecho de partículas só- 
lidas, que no suelen cambiar sus características. En este capítulo, se presenta el es- 
tudio del flujo de una sola fase fluida a través de una columna de partículas sólidas 
estacionarias. Un lecho de partículas se considera estacionario cuando todas sus 
características no varían. 


9.2. Ley de Darcy. Permeabilidad 


Después de una serie de experimentos Darcy demostró que la velocidad media, 
en una sección del lecho, era directamente proporcional a la caída de presión que 
experimenta el fluido al atravesar dicho lecho, e inversamente proporcional al es- 
pesor del mismo: 


L 


[ 9 . 1 ] 


en la que: (-A P) = Caída de presión a través del lecho. 

L = Espesor o altura del lecho. 

K = Constante de proporcionalidad. 

La constante K depende de las propiedades físicas del lecho y del fluido que 
circula a través del lecho. 

Esta ecuación indica que la relación entre la velocidad de circulación del fluido 
y la caída de presión que experimenta es lineal, lo que hace suponer que el régimen 
de circulación es laminar. Esto es así, ya que la velocidad a través de los intersti- 
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cios del lecho granular es baja, y además, la sección de cada uno de ellos es pe- 
queña; por lo que el valor del módulo de Reynolds no es elevado. 

Todo ello, presupone que la resistencia que ofrece el lecho al flujo del fluido es 
debida principalmente a rozamientos viscosos. Por lo que la ecuación anterior se 
suele expresar de la forma: 


1 (-A P) 

v = — - 

a r\L 


[ 9 . 2 ] 


en la que r\ es la viscosidad del fluido, mientras que la constante 1/a recibe el nom- 
bre de coeficiente de permeabilidad. 

La unidad de la permeabilidad es el Darcy , que se define como la permeabili- 
dad de un medio poroso al flujo viscoso para el paso de 1 ml/(s*cm 2 ) de un líquido 
con una viscosidad de 1 centiPoise, bajo una caída de presión de 1 atm/cm. 


9.3. Definiciones previas 

Los lechos porosos están constituidos por un conjunto de partículas sólidas, 
dispuestas generalmente al azar. Para la caracterización de la estructura de dichos 
lechos se deben definir una serie de variables, utilizadas en el estudio de la circula- 
ción de los fluidos a través de estos lechos. 

Superficie específica. También denominada área superficial específica, puede 
ser referida a todo el lecho o a la partícula. Para el lecho ( a s ) se define como el área 
de la superficie de lecho presentada al fluido por unidad de volumen del lecho. 

Área presentada al fluido . 

a s = expresada en m 1 

Volumen del lecho 

De modo análogo se puede definir el área superficial específica de las partícu- 
las (a s 0 ) como la relación entre el área de su superficie y su volumen. 

Área de partícula t 

a so = expresada en m 1 

Volumen de partícula 

En el caso de una partícula esférica de diámetro d n la superficie específica de 
partícula es a so = 6 ld r . 

Cuando las partículas no son esféricas, se define un diámetro equivalente de 
partícula (d P ), como aquel diámetro que poseería una esfera cuya relación área su- 
perficial a su volumen fuese la misma que la que posee la partícula. Este diámetro 
equivalente está relacionado con d r mediante un factor T, de tal forma que: 

dp — r d r 

en la que Y es el denominado factor de forma o esfericidad, cuyo valor depende 
de la forma de la partícula. Para el caso de partículas esféricas este factor de 
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forma es la unidad; sin embargo, para partículas diferentes este valor es menor 
que la unidad. En la tabla 9.1 se dan valores de la esfericidad para diferentes ti- 
pos de partículas. 

Además, el diámetro equivalente de la partícula se relaciona con su superficie 
específica según la expresión: 



Porosidad. El lecho de partículas no es compacto, sino que existen zonas li- 
bres de partículas, definiéndose la porosidad del lecho o fracción de huecos (e) 
como el volumen del lecho no ocupado por el material sólido. 

Debido a la porosidad del lecho, la superficie específica del lecho y de partí- 
cula no coinciden, sino que están relacionadas según la ecuación: 


fls = «5o(l-e) 


[93] 


En la bibliografía pueden encontrarse valores de a so y 8 para distintos tipos de 
lechos de partículas. Es fácil observar que cuando la porosidad aumenta, el flujo a 
través del lecho es mucho mejor, con lo que la permeabilidad aumenta. 

Para poder continuar estudiando las características del lecho, se va suponer que 
está compuesto por partículas dispuestas al azar, de forma que en el interior del le- 
cho se forman unos canales por los que circulará el fluido. Se supone que la longi- 
tud de estos canales (L’) es la misma para todos ellos, poseyendo el mismo diáme- 


Tabla 9.1. Esfericidad de partículas 


Forma de la partícula 


Esfericidad 


Esfera 

Cubo 


0,81 


Cilindros 
h = d 
h = 5d 


0,87 

0,70 

0,58 


h = 10 d 
Discos 


h = d/ 3 
h = d/6 
h = d! 10 


0,76 

0,60 

0,47 


Arena de playa 
Arena de río 
Distintos tipos de arena 
Sólidos triturados 
Partículas granulares 
Trigo 

Anillos Raschig 
Sillas Berl 


Tan alta como 0,86 
Tan baja como 0,53 


0,75 
0,5-0, 7 
0,7-0, 8 
0,85 

0,26-0,53 

0,30-0,37 


Fuente: Levenspiel (1993). 
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tro equivalente (. D e ). Todas las partículas están contenidas en una columna cilin- 
drica de diámetro D y cuya altura de lecho es L. 

A continuación se presentan una serie de variables, que se utilizarán para rela- 
cionar las diferentes características del lecho. 


Número de canales por m 2 de sección transversal de lecho: n’ 
Número total de canales del lecho: 

Área interfacial de 1 canal: 

Área interfacial del lecho: 

Volumen del lecho: 

Volumen de lecho ocupado por las partículas: V L (1 - e) 


n = n — D 2 
4 

n-D e U 

^-D 2 \ (nD e U) 


V l = — D 2 L 
L 4 


A partir de estas definiciones, es posible encontrar los valores de las superficies 
específicas: 


n’ — D 2 ) (7t D e U) «’ U n D e 

= \ 4 


L(l-e) 


[9.4] 


D 2 L( 1 - 8 ) 


\n , — D 2 )(TiD e U) n’ V n D 
gj= l 4 )_ 1 

— D 2 L 
4 


[9.5] 


por lo que es fácil correlacionar la superficie específica del lecho y de la partícula: 

a s = a so • ( 1 - e) 


Es posible definir el diámetro equivalente de un canal como cuatro veces el ra- 
dio hidráulico, siendo éste la relación entre el área de la sección de paso del canal 
y su perímetro mojado. 

Área de paso 

D e = 4R H = A —— F t 

Perímetro mojado 


D e = 4 


(Área de paso) Un 
(Perímetro mojado) Un 


4 (Volumen de un canal) n 
(Área interfacial de un canal) n 
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D 0 


A Volumen de huecos 

= 4 — 

Area interfacial total 


4 — D 2 Le 

4 

— D 2 La s 
4 s 


De esta forma, el diámetro equivalente de un canal será función de la fracción 
de huecos y de las superficies específicas, expresándose como: 


D = 


4 8 


4 8 


a s a so (1 ~~ e ) 


[9.6] 


Para futuras aplicaciones, es imprescindible correlacionar la velocidad de circu- 
lación del fluido a través de un canal (v c ) con la correspondiente al flujo a través de 
la columna libre de partículas (v). Para ello, se aplicará la ecuación de continuidad: 

p • S ■ V = p • S c ■ V C 

en la que p es la densidad del fluido, S y S c son las secciones de paso de la columna 
libre de partículas y para el conjunto total de canales, respectivamente. 

7t 

• Sección de paso de la columna: S = — D 2 

• Sección de paso de los canales: Será la suma de las secciones de paso de to- 
dos los canales, que como se han supuesto idénticos, se obtendrá multipli- 
cando el número de canales por la sección de paso de un canal. 


n 


S r = n — D: = n’ — D — D 


, K j^2 71 


2 

e 


Al introducir estas expresiones, de S y S c , en la ecuación de continuidad, es po- 
sible despejar la velocidad de circulación a través de un canal en función de la glo- 
bal, obteniendo la expresión: 

4 v 

v r = 

nn D e 

Pero si se tienen en cuenta las ecuaciones 9.4 y 9.6, es posible despejar de la 
9.4 el número de canales y de la 9.6 el diámetro equivalente del canal, que al susti- 
tuirlas en esta última expresión se obtiene: 

v c = v-f- [9.7] 

F 8 


expresión en que la velocidad de un canal es función de la global, y de las longitu- 
des del lecho y del canal, así como de la porosidad. 


248 Operaciones Unitarias en la Ingeniería de Alimentos 


9.4. Ecuaciones para el flujo a través de lechos porosos 

En el estudio de la circulación de fluidos a través de lechos porosos es de suma 
importancia el cálculo de la caída de presión que el fluido experimenta al atravesar 
dicho lecho de partículas. Esta pérdida de presión dependerá del régimen con el 
que circula el fluido, obteniéndose diferentes ecuaciones según sea el régimen la- 
minar o turbulento. 


9.4.1. Régimen laminar. Ecuación de Kozeny-Carman 

Si el flujo del fluido a través de los canales es laminar se podrá aplicar la ecua- 
ción de Fanning para cada uno de ellos: 

(-A P) _ 64 L’v 2 


p Re 2 D e 

pero si se tiene presente que el módulo de Reynolds es: 

P v c D e 


(Re) = 




y que la velocidad de circulación del fluido a través de un canal v c viene dada por 
la ecuación 9.7, se obtiene: 


(-A P) = 2 v r| 


( a s o) 2 (1 — £) 2 

C 3 


(L’f 

L 


se observa que la pérdida de presión depende, entre otros, de las longitudes de cada 
canal y del lecho. La longitud de cada canal U es superior a la del lecho. Si se su- 
pone que dichas longitudes son proporcionales U = K’L, y definiendo una cons- 
tante K” = 2 (X’) 2 , resulta: 


(-AP) = 2 v r| 


feo) 2 (1 ~ £ ) 2 

c 3 


(KTf 

L 


de la que se obtiene: 


(-AP) e 3 

tlL K”(l-e) 2 (a S0 ) 2 


[9.8] 


que se denomina ecuación de Kozeny-Carman. 

Al comparar esta ecuación con la de Darcy (ecuación 9.2), se obtiene que la 
permeabilidad será: 

1 1 e 3 


a K” (1 - e) 2 (a S0 ) 2 
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La constante K ” se denomina constante de Kozeny. En lechos en los que la po- 
rosidad y superficie específica no varían con el espesor del lecho, se ha encontrado 
experimentalmente que esta constante posee un valor de 5 ± 0,5. En realidad su va- 
lor depende del tipo de relleno y de la porosidad, tomando distintos valores según 
sea la forma de las partículas y porosidad del lecho. En el caso que las partículas 
sean de forma esférica, el valor de esta constante es de 4,8 ± 0,3. En la bibliografía 
se pueden encontrar valores de esta constante para diferentes tipos de relleno. 

El valor de la constante de Kozeny no es el mismo para todos los tipos de re- 
lleno, sino que depende de la relación L’/L. Carman ha demostrado que: 


K” = K 0 (L7L) 


La relación L’/L recibe el nombre de tortuosidad, mientras que K 0 es un factor 
que depende de la sección transversal del canal. A pesar de que la tortuosidad y el 
factor K 0 pueden variar, esta variación es tal que cuando una aumenta la otra dis- 
minuye, y viceversa, de forma que su producto posee valores próximos a 5. 

En la práctica se observa que en los lechos de partículas, las que se hallan en 
contacto con la pared están menos compactadas, lo que implica que la resistencia 
que ofrece el lecho al flujo es menor que el dado por la ecuación de Kozeny Car- 
man. De forma experimental Coulson ha obtenido un factor de corrección K P , de 
modo que se tenga presente este efecto. Dicho valor viene dado por la ecuación: 



[9.9] 


en la que A P es la superficie de la pared de la columna que contiene el lecho por 
unidad de volumen de dicho lecho 

Para el cálculo de la caída real de presión, deberá multiplicarse la calculada a 
partir de la ecuación de Kozeny-Carman por el valor de este factor K P \ 


(-AP) real = K p (-AP) 


9.4.2. Régimen turbulento. Ecuación de Burke-Plummer 

En el apartado anterior se ha supuesto que el fluido circula por los canales de 
forma laminar, pero la ecuación obtenida no es válida en muchos casos. Por ello, si 
el régimen de circulación es turbulento puede pensarse en aplicar la ecuación de 
Fanning a la circulación del fluido por un canal; así: 


(-A P) (v c ) 2 

P 7 2 D e 


teniendo en cuenta las expresiones del diámetro equivalente y de la velocidad por 
un canal [ecuaciones 9.6 y 9.7]: 


(-A F) 1 (v¿’) 2 ¿'flaQ-e) 

p 7 2 (eL) 2 4 £ 
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Al igual que antes, si se supone que U = K’L, y que la superficie específica está 
relacionada con el diámetro de la misma por la expresión a so = 6 !d P , se obtiene: 


(-A P) 
P 


= 3 (K’ff 


p V 2 (1 -8) 

d p £ 3 


si se define un factor de fricción modificado/’ =f(K’) 3 se obtiene la ecuación: 


(-A P) p V 2 (1 — 8) 

p 1 d p 8 3 


[9.10] 


Esta ecuación se denomina de Burke-Plummer, en la que el valor de/’ se ob- 
tiene a partir de experimentación, dependiendo del número de Reynolds. En el 
apartado siguiente se verá el modo de obtener el valor del factor de fricción modi- 
ficado. 


9.4.3. Flujo global laminar-turbulento. Ecuaciones de Ergun 
y Chilton-Colburn 

Hasta ahora se han obtenido ecuaciones para la circulación en régimen lami- 
nar y turbulento, por separado. Por tanto, sería conveniente obtener una sola 
ecuación que se pudiera utilizar para los dos tipos de regímenes de circulación a 
la vez. Observando las ecuaciones de Kozeny-Carman y de Burke-Plummer, 
puede pensarse que la pérdida de presión por unidad de longitud del lecho será 
una expresión del tipo: 


L 


o bien una combinación lineal de dichas ecuaciones: 


(-A P) _ , 36 K” (1 - 8) 3 
L ~ a eUp 


. Q , 3/’(l-e)p 2 
V + P ~ eX V 


[9.11] 


Para la obtención de las distintas constantes (a’ y (3’) que aparecen en esta 
ecuación se recurre a datos obtenidos por experimentación, ajustándolos a ella. 

En la circulación de fluidos por el exterior de sólidos, no sólo se utiliza el mó- 
dulo de Reynolds para un canal, sino que en algunos casos se suele modificar, o 
bien se utiliza el módulo de Reynolds de partícula. A continuación se dan las ex- 
presiones para cada uno de ellos: 

_ , , . , _ p v c D e 4p vU 

• Reynolds para 1 canal: Re = = 

T| a so (1 - 8) T) L 


pv 

aso (1 - e) r| 


Reynolds modificado: Re’ = 
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• Reynolds de partícula: 


Re n 


P vdp 

r| 


Cuando el valor del módulo de Reynolds de partícula es menor de 40 (Re p < 40) 
el régimen de circulación se considera laminar, mientras que para valores superio- 
res (Re p > 40) el régimen se considera turbulento. 

De la ecuación 9.10, es fácil obtener la expresión para el factor de fricción mo- 
dificado: 


(-A P)d p e 3 

L (1 - 8) p V 2 


[9.12] 


Este factor de fricción suele obtenerse a partir de gráficas, en las que 3/’ se re- 
presenta frente a Re p /(1 - 8), en coordenadas doble logarítmicas (figura 9.1). 

También puede obtenerse en otro tipo de gráfica, en la que/72 se representa 
frente al módulo de Reynolds modificado (Re’), también en coordenadas doble lo- 
garítmicas (figura 9.2). 

A partir de una serie de datos experimentales, de valores del factor de fricción 
modificado 3/’, calculados a partir de la ecuación 9.10, se obtienen los valores co- 
rrespondientes de sus Reynolds de partícula. Se busca una función que correla- 
cione 3 f’ con Re p : 

3/'=<|>(Re„) 


Del ajuste de los datos experimentales la ecuación obtenida es la siguiente: 


3 /’ = 


150 (1-e) 
Re ;) 


+ 1,75 



Figura 9.1. Factor de fricción modificado (3f’) en función del módulo de Reynolds 
de partícula (Re p ) ( Adaptado de Foust et al., 1980). 
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Figura 9.2. Factor de fricción modificado (f’) en función del módulo de Reynolds modificado 
(Re’) ( Adaptado de Coulson y Richardson, 1981). 


Si en esta ecuación se sustituyen los valores de 3/\ dados por la ecuación 9.12, 
y el correspondiente al módulo de Reynolds de partícula, se obtiene la expresión: 


(-AP) 

L 


= 150 


(1 - e) 2 r| 

£ 3 d,; 


v + 1,75 


(1 — e) p 

£ 3 d p 


[9.13] 


Esta última expresión se denomina ecuación de Ergun, y puede utilizarse para 
el cálculo de la pérdida de presión que experimenta el fluido al atravesar un lecho 
relleno, independientemente de cual sea el tipo de régimen de circulación. 

Al comparar las ecuaciones 9.11 y 9.13, fácilmente se observa que: 

36 K” a’ = 150 


3/’ P' =1,75 


Cuando en una columna rellena de partículas circulan en contracorriente un gas 
con un líquido, es conveniente utilizar la denominada ecuación de Chilton-Col- 
burn, que es una ecuación empírica basada en la de Fanning: 

2f'^~ [9.14] 

L a p 

El factor de fricción modificado/’ puede calcularse a partir de las figuras 9.1 ó 
9.2, o bien pueden utilizarse las siguientes expresiones: 

• Para régimen laminar (Re p <40): /’ = 850/Re p 

38 

• Para régimen turbulento (Re p > 40): /’ = 

(R e p ) ’ 
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La ecuación de Chilton-Colburn se puede utilizar cuando el lecho está formado 
por partículas macizas, pero cuando las partículas son huecas, debe multiplicarse el 
segundo miembro de la ecuación por un factor K r , que viene dado por la ecuación: 

K _ 0.24 
r ( d P ) m 

en la que d P es el diámetro nominal de las partículas expresado en pulgadas. 

Los huecos creados en las inmediaciones de la pared también pueden afectar al 
cálculo de la caída de presión, siendo necesario introducir un nuevo factor de co- 
rrección de pared K P , que depende del régimen de circulación y de la relación diá- 
metro partícula a diámetro del lecho (d P /D). En el caso que esta relación sea menor 
de 1/6 no se considera el efecto de pared. Para el cálculo de este nuevo factor K P se 
utiliza la figura 9.3, en la que K P se halla representado frente a la relación d P /D\ 
existiendo dos curvas, una para régimen laminar y otra para turbulento. 

Por tanto, en general, la caída real que experimentará el fluido se obtendrá mul- 
tiplicando la caída de presión calculada a partir de la ecuación de Chilton-Colburn 
por estos dos factores; es decir: 


( ~^P)Real ~ K r K p (“A P) Teórico 



Figura 9.3. Factor de corrección del efecto de pared (K P ) en función de la relación diámetro 
partícula/diámetro de columna (d/D) (adaptado de Sawistowski y Smith, 1967). 
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9.5. Fluidización 

Cuando un fluido circula a través de un lecho de partículas, si su velocidad de 
circulación es baja, el lecho permanece estático. Pero si su velocidad aumenta, 
puede ocurrir que el lecho se expanda, existiendo una reordenación de las partícu- 
las, con el consiguiente aumento de la porosidad del lecho. 

Si se supone un lecho de partículas, a través de las cuales circula un fluido al 
que se le va aumentando paulatinamente la velocidad de circulación. A velocidades 
bajas el lecho permanece estático, sin aumento de su altura, por lo que al represen- 
tar en coordenadas doble logarítmicas la caída de presión frente a la velocidad de 
circulación se obtiene una recta, tal como la representada en la figura 9.4 por el 
segmento OA. 

Si se aumenta más la velocidad, las partículas empiezan a separarse, pero toda- 
vía permanecen en contacto, aunque la relación entre la caída de presión y la velo- 
cidad continúa siendo lineal, pero la pendiente es menor (segmento AB). En el 
punto B, en el que las partículas no se hallan en contacto se dice que el lecho es 
fluidizado. A partir de este punto, si se aumenta la velocidad puede haber una pe- 
queña caída de presión (segmento BC), pero si se sigue aumentando la velocidad, 
la caída de presión aumenta linealmente con la velocidad, pero con una pendiente 
menor, hasta que la velocidad es lo suficientemente elevada como para arrastrar las 
partículas, cosa que ocurre a partir del punto D. 

La velocidad para la cual ocurre la fluidización se conoce como mínima de 
fluidización; mientras que la correspondiente al punto de arrastre es la velocidad de 
arrastre. 

En el tramo en que el lecho está fluidizado, la velocidad del fluido en contacto 
con las partículas es superior a cuando las ha superado, ello hace que las partículas 



Figura 9.4. Caída de presión en lechos fluidizado s. 
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no queden arrastradas, sino que vuelven a caer el lecho. El conjunto de partículas 
posee un movimiento desordenado, presentando un aspecto como si las partículas 
estuvieran en ebullición, denominándose lechos hirvientes o fluidización disconti- 
nua a este tipo de proceso. Cuando las partículas son arrastradas, el tipo de fluidi- 
zación se denomina continua, y es un tipo de circulación de dos fases, que es la 
base del transporte neumático. 


9.5.1 . Velocidad mínima de fluidización 

La velocidad a la que empieza la fluidización del lecho se conoce como veloci- 
dad mínima de fluidización. En este punto existe un equilibrio dinámico entre la 
fuerza que el campo gravitatorio y el fluido ejercen sobre las partículas. 

La fuerza gravitatoria ejercida sobre las partículas viene dada por la expresión: 

F G = (p P -p)SL(l-e)g [9.15] 

en la que p P y p son las densidades de las partículas y del fluido, respectivamente; 
S la sección de paso de la columna que contiene las partículas; L la altura del lecho; 
8 su porosidad, y g la constante de gravedad. 

La fuerza que el fluido ejerce sobre el lecho de partículas es la de presión, y que 
se calcula por la expresión: 

F P = (-AP)S [9.16] 


en la que (-A P) es la caída de presión que experimenta el fluido al atravesar el le- 
cho, y cuya expresión depende del régimen de circulación del fluido. 

Por tanto, para el cálculo de la velocidad mínima de fluidización deben igua- 
larse las dos fuerzas, la de gravedad y la de presión F G = F P , o lo que es lo mismo: 

(p P -p)SL(l-e)g = (-AP)S [9.17] 


La expresión de la caída de presión viene dada por la ecuación de Ergun: 


(-A P) 
L 


= 150 


(1 -8) 2 T| 

e 3 # 


v+ 1,75 


(1 - e) p 

£ 3 d p 


V 


2 


A pesar de que ésta es la ecuación general para el cálculo de la caída de presión, 
según sea el régimen de circulación del fluido podrá simplificarse. A continuación 
se obtendrán las diferentes ecuaciones que permitirán el cálculo de la velocidad 
mínima de fluidización, dependiendo del tipo de régimen de circulación del fluido. 


9.5 . 1. 1. Régimen laminar 

Cuando el régimen de circulación del fluido sea laminar, el segundo término 
del segundo miembro de la ecuación de Ergun es despreciable frente al primero, 
con lo que su expresión se simplifica: 
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(-A P) 

L 


= 150 


(1 -e) 2 r| 
e 3 dp 


v 


por lo que al sustituir esta expresión en la ecuación 9.17, es posible determinar la 
velocidad mínima de fluidización en régimen laminar 


V mf 


1 


(£ m /) 3 


Pp~P 


150 (l-8 m/ ) T) 


gd, 


[9.18] 


9.5. 1.2. Régimen turbulento 


Cuando el régimen de circulación es turbulento, en la ecuación de Ergun, es el 
término de velocidad el que es despreciable frente al del cuadrado de la velocidad. 
Por tanto, la caída de presión que experimenta el fluido se podrá expresar como: 


(-AP) 

L 


1,75 


(l-e)p 
£ 3 d p 


V 


2 


Al sustituir esta expresión en la ecuación 9.17, es posible determinar la veloci- 
dad mínima de fluidización en régimen turbulento: 


V mf 


= 0,756 


P,-P 


g (£m/) 3 d p 


1/2 


[9.19] 


9.5 . 1.3. Régimen de transición 


Puede ocurrir que el fluido circule con un régimen superior al laminar, pero no 
esté completamente desarrollado el turbulento. En estos casos es necesario aplicar 
la ecuación de Ergun para el cálculo de la caída de presión. A pesar de ello, se 
aconseja utilizar esta ecuación cualquiera que sea el régimen de circulación. 

En estos casos, al sustituir la expresión de Ergun en la ecuación 9.17 se obtiene 
una ecuación de segundo grado, que es necesario resolver para hallar la velocidad 
mínima de fluidización: 


(1 - e m/ ) (p p - p) g = 150 


( £ m/) dp 


v mf + 1,75 


(1 ~ £ m/) P 

( £ m/) dp 


(v m/ ) 2 


[9.20] 


9.5.2. Porosidad mínima de fluidización 

Lo mismo que ocurría con la caída de presión del fluido, que variaba con la ve- 
locidad de circulación, la porosidad del lecho también experimenta variación. Para 
velocidades bajas, la fuerza de presión no es suficiente para variar la estructura del 
lecho, y su porosidad no varía, pero a velocidades más altas, el lecho se expansiona 
aumentando el volumen de huecos. 
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Si se representa en coordenadas doble logarítmicas la variación de la porosidad 
con la velocidad lineal de circulación se obtiene una gráfica como la representada 
en la figura 9.5. A velocidades bajas la porosidad no varía, pero llega un momento 
en que aumenta con la velocidad, siendo esta variación lineal. La porosidad en el 
momento que empieza la fluidización no se corresponde al punto en que la porosi- 
dad empieza a aumentar, sino que ocurre a una velocidad más alta, que es la mí- 
nima de fluidización. 

Para el cálculo de la porosidad mínima de fluidización pueden utilizarse unas 
gráficas, que existen para cierto número de materiales (McCabe y Smith, 1968). 
Sin embargo, cuando no se dispone de datos, es posible utilizar la expresión: 


e mf = 1 - 0,356 (log d P - 1) [9.21] 

ecuación en la que el diámetro de partícula debe expresarse en micrones (10 -6 m), 
y es válida para valores de d P comprendidos entre 50 y 500 micrones. 

9.5.3. Altura del lecho 

Evidentemente, al aumentar la velocidad de circulación del fluido a través del 
lecho poroso, no sólo aumenta la caída de presión, sino que el aumento de la poro- 
sidad es debido a que la altura del lecho también aumenta. Este aumento de altura 
de lecho está íntimamente ligado al aumento de porosidad, de forma que las poro- 
sidades correspondientes a dos alturas cualesquiera están relacionadas por la ex- 
presión: 

Ml-£i) = L 2 (l-e 2 ) [9.22] 



Figura 9.5. Porosidad de lechos fluidizados. 
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De forma particular, si para el lecho fijo le corresponde una porosidad 8 0 y una 
altura de lecho L 0 , la altura y porosidad para otro instante están relacionadas con 
ésta por la ecuación: 

L(l-£) = L 0 (l-£ 0 ) 

Si el lecho fuese compacto, las partículas ocuparían todo el lecho y no existi- 
rían huecos, por lo que la expresión que correlaciona las alturas se simplificaría: 

L(1-e) = L c 

en la que L c es la altura del lecho compacto. 


